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ABSTRACT 


A  simple,  approximate  method  is  presented  for  rapid  determination 
(utilizing  mini-computers)  of  plume  shapes  produced  by  jet  expansion 
fr'-’i  axially  symmetric  nozzles.  The  analysis  is  based  on  concepts 
developed  by  Johannesen  and  Meyer.  It  is  shown  how  the  method  for  plume 
shape  determination  can  be  utilized  to  model  at  least  the  geometric 
aspects  of  prototype  plumes  as  well  as  to  accour.t  for  significant 
inviscid  and  viscid  aspects  of  the  base  flow  problem.  A  calculation 
procedure  and  numerical  examples  are  presented  to  guide  the  reader  in 
applying  the  modeling  technique. 
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SYMBOLS 


a* (ft/sec) 

critical  acoustic  velocity 

c1  (-) 

constant  of  integration  in  Eq.  (5) ,  determined  from 
Eq.  (13) 

M(-) 

Mach  number 

M*  (-) 

critical  Mach  number 

R(-) 

polar  coordinate  radius  vector  (dimensionless)  (exit 
radius  of  nozzle  =  1) 

rc("> 

initial  radius  of  curvature  of  expanded  free  jet 
boundary 

u(£ c/sec) , 
v(ft/sec) 

velocity  components 

O 

O 

constant  of  integration  in  Eq.  (6) ,  determined  by 

Eq.  (14) 

Y  (“) 

specific  heat  ratio 

<M°) 

polar  angle 

A(-) 

[  (Y  ~  D/(Y  +  1)]% 

n(-) 

A  (<f>  +  Bq)  ,  auxiliary  angle 

e  (°) 

streamline  angle 

u(°> 

Mach  angle 

w(#) 

Prandtl-Meyer  streamline  turning  angle  subscript 

Subscripts 

F 

conditions  at  final  expansion  fan  line  as  R  +  0 

L 

conditions  at  initial  expansion  fan  line  as  R  +  0 

3 

free  jet  boundary  as  R  +  0 

M 

model 

P 

prototype 

iv 


I.  INTRODUCTION 


Plume-slipstream  interactions  are  of  importance  in  propulsion, 
stability,  and  guidance  problems  related  to  missile  flight.  Depending 
on  design  performance  and  flight  profile,  low  pressures  in  base  regions 
may  impair  propulsive  efficiency  by  imposing  drag  penalties  (base  drag) 
and  also  generate  there  an  unfavorable  thermal  environment  (base  heating) . 
With  increasing  jet-to-ambient  pressure  ratio,  the  plume  size  increases, 
adverse  base  drag  effects  gradually  disappear,  and  larger-than-ambient 
base  pressures  produce  base  thrust.  This  can  eventually  lead  to  flow 
separation  from  the  missile  afterbody,  a  situation  which  may  be 
considered  beneficial  in  view  of  possible  drag  reduction  but  disturbing 
if  bnffetting  and  loss  of  control  result  from  such  interactions.  The 
combined  effects  of  plume  induced  separation  and  angle  of  attack  of 
the  missive  can  give  rise  to  destabilizing  moments  [l]. 

Although  comprehensive  computer  programs  allow  the  analysis  of  jet- 
slipstream  interactions  for  supersonic  flight  velocities,  they  fail  to 
give  coverage  for  the  transonic  flight  regime  which  is  most  vulnerable 
to  the  adverse  effects  of  plume  induced  separation.  Consequently,  the 
results  of  well  planned  experiments  continue  to  serve  as  the  main  source 
for  guidance. 

The  observation  that  transonic  flows  past  afterbodies  are  not 
capable  of  large  angular  deflections  in  negotiating  adverse  pressure 
gradients  places  much  emphasis  on  the  study  of  plume  boundaries  as  they  > 
establish  geometrical  constraints  on  the  transonic  flow  approaching  the 
base. 


The  analysis  of  plume  shapes  as  they  are  affected  by  nozzle 
configuration,  propellants,  and  pressure  ratio  attracts  much  interest. 
Utilization  of  well  established  computer  programs  based  on  the  method 
of  characteristics  for  axially  symmetric  supersonic  flows  is  here 
supplemented  by  a  simple  MINI-COMPUTER  oriented  analysis  based  on 
earlier  work  by  Johannesen  and  Meyer  [2]. 


1 


II.  OBJECTIVES 


The  method  presented  here  shall  allow  the  rapid  determination  of 
plume  shapes  produced  by  jet  expansion  from  axially  symmetric  nozzles. 
The  effects  of  nozzle  geometry  (area  ratio--or  nominal  exit  Mach  number 
and  nozzle  wall  shape — divergence  angle)  propellant  composition  ( y , 
specific  heat  ratio)  and  overall  pressure  ratio  (M^,  or  (p/P  )  |  )  are 

to  be  considered.  F 

The  approximate  analytical  method  shall  furnish  a  rational  basis 
for  plume-modeling. 
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III.  ANALYSIS 


The  analysis  follows  the  concepts  developed  by  Johannesen  and 
Meyer  [2]  in  expressing  the  flow  field  near  the  centered  expansion  at 
the  nozzle  exit  in  the  form  of  series  expansions  with  respect  to  the 
radius  vector  R  (Fig.  1).  The  velocity  components  thus  are: 

u  =  u  (<p)  +  u.  (4>)  R  +  CHR2)  (1) 

o  1 

v  =  v  (<}>)  +  v  (4>)  R  +  0(R2)  (2) 


and  after  substitution  into  the  conservation  equations  for  axially 
symmetrical  potential  flow,  one  arrives  at  a  system  of  equations  for  the 
yet  unknown  functions  uq($),  u^($),  vq(<J>),  and  v^(<J>)  which  can,  after 

comparing  terms  of  equal  order,  be  solved  for  given  initial  and  final 
boundary  conditions.  Two  types  of  solutions  appear  for  the  two  types 
of  flow  regions  A  and  B  (Fig.  2) : 

1)  For  the  centered  expansion  <fc  <  $  <  <j>  (Region  A  is  identified 

L  r 

in  Fig.  2,  a  modified  Prandtl-Meyer  fan) ,  one  obtains 


u  (4>)  ■, 

_?  =  £  sin  lxw  +  6  )  ] 

d  A  O 

v  <♦> 

- J-  =  cos  [ X (<f>  +  B  )  ] 


(4>) 


i(cos  n)^3Y-l]/[2(Y-l)]}  (sin  n)% 


(3) 

(4) 


•  j[l1  (X)  -  Xl3(n)  ]  cos  Bo 


+  [i2  (n)  +  u4(n)  ]  Sin  6o  +  c1 


where 


n  =  X(<J>  +  6  ) 
o 


and 


Ix(n) 


/ 


cos  (  (sin  n)  2  (c 


-  (1/2X2)  . 
3D)  dn 


(5) 


(6) 
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i2  (n) 


1 1 

f 


sin(Y)(sin  n)  (cos  n) 


-(1/2X2) 


dn 


Vn) 


=  j  sin  (?)  (si 


,  •  .  -3/2  .  .  { [y-3]/[2  (y-1)  ]  } 

(sin  n)  (cos  n)  dn 


Vn) 


■  /  °“  (») 


(sin  n)'3/2  (cos  n)  t>-3]/[2(v-l)]  an 


vL  (<f>> 


1  \ 

2  a* 


(7) 


while 


,  4  [!a  sin  ♦  -  ^1  +iiC-23LzJi(^5|2 

U1  a  Lf  a*  11  a  a  L  Y  +  1  \V0/J 


u 

2  -°~ 


(8) 


2)  Outside  the  centered  wave  region  one  arrives  at  relations  which 

allow  establishing  boundary  conditions  from  matching  solutions  at  the 

lower  (<}>  )  and  final  (4>_)  Mach  lines  of  the  centered  fan. 

L  r 

(a)  The  approaching  flow  yields  (Region  B  is  identified  in 
Fig.  2  and  shown  in  Fig.  3) . 


u  (<L) 

O  L 


=  M*  cos  u 

L  Jj 


(9) 


v  (<i>  ) 

O  L  . 

sln  \ 


VV  ( 

r*  =  mt.  sin  2y- 


(10) 


Toe 

‘L  jj3R 


_1_  dM^ 

_  .  M*  dr 
R=0  L 


cot 


R=0 


-1 


sin  0  sin2  u 
L  L 


(id 


accounting  for  the  flow  conditions  upstream  of  the  nozzle  lip,  in 
particular,  for  the  wall  curvature  (d0/dR) I and  the  flow  acceleration, 
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R-*-0 


Figure  3.  Boundary  Conditions  Imposed  by  Approaching  Flow 

For  a  conical  nozzle,  where  dG/dR  =  0,  the  assumption  of  source  flow 
yields 


u  (<|>  )  4  sin  0  cos2  t 

-L~L-  =  M*  -  L  L 


,  >2  S8 

Y  “  1  1  -  X2  M* 
L 


(12) 


so  that  in  Eq.  (5)  can  be  evaluated 


C,  = 


8  XM*  sin  0  cos2  u 

L  Li  Li 


-  (13) 


X2  M*2 

Li 


Y  "  1  1  -  X2  M*2 
L 


(cos  n  ) 
Li 


{ [ 3y-l ]/[ 2  (y-1)  ] } 


(sin  nL) 


The  constant  8  us  determined  throuqh 
o 


i  .  -i  x 

8  =  -  tan  - - — 

o  X  tan 


(14) 
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Figure  4.  Boundary  Conditions  Imposed  by  Expanded  Flow 


where 


-  G  -  u  +  90°; 


G  =  9  +  (D  (M  )  -  w  (M  ) 

r  U  £  Li 


(17) 


For  the  initial  curvature  of  a  free  jet  (plume)  where  <3M*/dR  =  0, 
one  obtains,  therefore. 


d0 

dR 


R+0 


-v  r- 

sin  2 ^ 


sin  9_  sin2  y  + 

r  F 


uiV.  _i_” 

a*  M* 
F_ 


(18) 


The  initial  of  curvature  is 


1 


ae 

dR 

R-H3 

Obviously,  Eq.  (16)  can  also  be  used  to  determine  the  initial  Mach 
number  (or  pressure)  distribution  along  a  v?all  boundary  to  which  the 
nozzle  flow  expands.  In  this  case  (d0/dR)  |  and  0p  are  considered  to 

be  given  and  dM*/dR  is  to  be  solved  for. 
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IV.  PLUME  MODELING 


For  any  given  input  and  Mp) ,  there  results  an  approximate 

plume  shape  defined  by  the  initial  slope  of  the  jet  boundary  0  and  the 

(dimensionless)  radius  of  curvature  r  .  F 

c 

1.  MODELING  LAWS 

a.  Geometric  Modeling 

If  plume  modeling  is  to  be  achieved,  it  will  be  necessary  to 
match  geometrically  (within  our  degree  of  approximation) 


and 


(20) 


One  expects  also  that  the  following  information  is  given  for  the 
prototype : 


VTop'W9Lp'  and  MFp  • 

For  the  model  it  shall  be  assumed  that  the  propellant  gas  has  been 
selected  so  that  Rm  and  ym  are  determined.  As  we  consider  0  ,  M  ,  and 

lm  lm 

M  as  the  parameters  to  be  determined,  the  geometrical  constraints 
M 

(Eqs.  (19)  and  (20))  must  be  supplemented  by  one  additional  condition: 
The  proper  modeling  of  the  closure  conditions  for  the  wake  is  necessary. 

b.  Specifying  Model  Laws 

Ideally,  this  specifying  condition  should  properly  account  for 
the  viscid  aspects  of  the  base  flow  problem  in  their  interaction  with 
the  inviscid  components.  One  must,  on  the  other  hand,  realize  that  the 
comprehensive  nature  of  the  componenc  model  for  jet-slipstream  inter¬ 
actions  in  the  vicinity  of  propulsive  afterbodies  [3]  which  apparently 
reflects  physical  reality,  in  essence  requires  two  closure  conditions 
to  be  satisfied.  The  first  relates  to  the  recompression  ratio  at  the 
end  of  the  wake  (a  dynamic  condition) ,  and  the  second  condition  calls 
for  conservation  of  mass  in  the  wake.  With  only  one  choice  to  be  made 
for  the  closure  of  the  model  law,  one  must  raise  the  question  whether  or 
not  a  logical  choice  can  be  made  so  that  modeling  can  be  accomplished 
in  the  most  constructive  manner.  If  this  should  prove  impractical,  one 
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could  address  oneself  again  to  the  comprehensive  analytical  model  to 
augment  the  experimental  modeling  procedure,  e.g. ,  by  a  theoretically 
predetermined  base  bleed. 


With  these  reservations,  we  shall  now  consider  a  number  of  simple 
specifying  conditions  for  closure  of  the  modeling  law: 


(1)  If  one  adopts  for  this  purpose  the  simple  recompression 
model  of  the  Chapman-Korst  model  (expecting  that 
empirical  corrections  cancel  out  between  model  and  proto¬ 
type  flows,  it  follows  that 


(X  -  *4  t)  <V/ V1>  '  i1  ' 


A  M*‘ 


2  \(Yj/(Y  -1) 


d  "p  F 
P  P 


which,  fo  *  (}i2  AM*2  <<  1  and  ^  reduces  to 

aM  p 


M*  =  M* 
F„  F 
M  p 


'M 


Yp  +  V  YM  +  1 


(21a) 


(2)  Matching  the  momentum  pu2  of  corresponding  plume  boundaries 
yields,  since  the  pressure  along  plume  surfaces  are  to  be 
equal , 


(3)  Matching  the  flux  density  pu  would  result  in  the  implicit 
relation  - - 


M* 


M* 


M 


1  -  A2  M*2  1  -  V2  M* 


M  F 


M 


p  FpN 


y 


(Y  +  1)  R  T 


E. 


(Y  +  1)  R  T 


M 


which  can  be  solved  as  a  quadratic  equation  in  M*  .  It 

FM 

is  interesting  to  note  that  there  exists  the  possibility 
to  explore  the  effect  of  stagnation  temperature  variation, 

(4)  Matching  the  supersonic  inviscid  streamline  deflection- 

pressure  rise  relation  on  the  basis  of  local  linearization 
(weak  shock  approximation)  requires 


Y 


P 


V 


1 


Y  M2 
rM  M 


(21d) 
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This,  like  other  approximations,  will  impose  certain 
mathematical  and  physical  limits  on  the  range  of 
possible  modeling. 

2.  CALCULATION  PROCEDURE 

We  intend  to  model  a  prototype  plume  which  results  from 


V  W  *”d  % 


yielding  an  initial  slope  0_  and  an  initial  radius  of  curvature  r 

F  cp 

For  a  given  model  gas  (R  ,T  „,y..)  selection  of  the  proper  specifying 

M  OM  M 

condition,  e.g.,  Eq.  (21a),  will  determine  M*  . 

F 

M 


Combining  Eqs.  (17)  and  (19)  yields 

0..  =  6T  +  u(M_  » Y  )  -  u(M_  ,Y  )  -  m(M_  ,ym)  +  w  (M  ,ym)  (22) 

L,,  L_  F  p  L  p  F  M  L„  M 

MP  P  p  M  M 

where  0  and  M  are  unknown. 

Xj  Xj 

M  M 


For  any  selected  trial  value  of  (there  results  now  a  0^ 


that  the  calculation  procedure  outlined  in  Section  3  will  yield  the 
curvature  r  .  If  a  solution  exists,  it  will  be  found  by  satisfying 

'■'w 


r  =  r  (Fig.  5).  It  is  noteworthy  that  nozzle  geometries,  nozzle 
CM  S 


Mach  numbers,  and  jet- to- ambient  pressure  ratios  can  differ  widely  between 
model  and  prototype,  particularly  also  as  a  consequence  of  changing  the 
specifying  condition.  On  the  other,  even  these  wide  parametric  variations 
do  not  seem  to  impair  the  accuracy  of  matching  the  plume  surface  geometry. 
Thus,  the  present  method  affords  a  convenient  approach  to  identifying 
and  establishing  the  usefulness  of  the  best  suited  specifying  condition 
so  that  the  entire  plume-slipstream  interaction  is  properly  modeled. 


3.  NUMERICAL  EXAMPLE 


/ 


The  prototype  plume  is  defined  by 


Y 


P 


1.2 
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2 


0.02126  . 


Modeling  is  to  be  done  by  air 
%  -  l-4  • 

Selection  of  different  specifying  conditions  produces  the  following 
parameters : 


TABLE  I 


MODEL,  Eq.  (21a) 

MODEL,  Eq.  (21b) 

Y 

1.2 

1.4 

1.4 

ml 

2.00 

3.515 

1.394 

0L 

6° 

23.6° 

1.75° 

mf 

3.00 

4.651 

2.7775 

6F 

38.195° 

38.2° 

38.2° 

r 

4.2788 

4.27 

4.28 

c 

Shown  in  Figure  6  are  the  plume  contours  for  the  prototype  flow 
and  the  two  modeled  flows  as  produced  by  the  method  of  characteristics 
[3].  Agreement  is  excellent  over  the  whole  calculated  range  which 
assures  the  usefulness  of  the  geometrical  aspects  of  the  modeling 
procedure  independent  of  the  selection  of  a  specifying  condition.  Also 
plotted  in  Figure  6  is  the  circular  arc  approximation  which  forms  the 
basis  for  the  modeling  procedure.  The  circular  arc  approximation,  in 
its  own  rights,  produces  reasonable  plume  boundaries  within  the  range 
of  one  nozzle  exit  radius. 


13 


REFERENCES 


t 


1.  Deep,  R.  A.,  J.  H.  Henderson,  and  C.  E.  Brazzel,  "Thrust  Effects  on 
Missile  Aerodynamics,"  US  Army  Missile  Command,  Redstone  Arsenal, 
Alabama,  Report  No.  RD-TR-71-9,  May  1971. 


Johannesen,  N.  H.  and  R.  E.  Meyer,  "Axially-Symmetrical  Supersonic 
Flow  Near  the  Centre  of  an  Expansion,"  The  Aeronautical  Quarterly, 
Vol.  2  (1950),  op.  127-142. 


3.  Addy,  A.  L. ,  "Analysis  of  the  Axisymmetric  Base-Pressure  and  Base- 
Temperature  Problem  with  Supersonic  Interacting  Freestream-Nozzle 
Flows  Based  on  the  Flow  Model  of  Korst,  et  al.,  Part  III:  A 
Computer  Program  and  Representative  Results  for  Cylindrical, 
Boattailed,  or  Flared  Afterbodies,"  US  Army  Missile  Command, 
Redstone  Arsenal,  Alabama,  Report  No.  RD-TR-69-14,  February  1970. 


15 


